Effects of Dirac sea on pion propagation in asymmetric nuclear matter 
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We study pion propagation in asymmetric nuclear matter (ANM). One of the interesting con- 
sequences of pion propagation in ANM is the mode sphtting for the different charged states of 
pions. First we describe the pion-nucleon dynamics using the non-chiral model where one starts 
with pseudoscalar (PS) ttN coupling and the pseudovector (PV) representation is obtained via suit- 
able non-linear field transformations. For both of these cases the effect of the Dirac sea is estimated. 
Subsequently, we present results using the chiral effective Lagrangian where the short-distance be- 
havior (Dirac vacuum) is included by re-defining the field parameters as done in the modern effective 
field theory approach developed recently. The results are compared with the previous calculations 
for the case of symmetric nuclear matter (SNM). Closed form analytical results are presented for the 
effective pion masses and dispersion relations by making hard nucleon loop (HNL) approximation 
and suitable density expansion. 
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I. INTRODUCTION 

Pions in nuclear physics assume a special status. It 
is responsible for the spin-isospin dependent long range 
part of the nuclear force. In addition, there are variety 
of physical phenomena related to the pion propagation 
in nuclear matter. One of the fascinating ideas in 
relation to the pion-nucleon dynamics in nuclear matter 
is the pion condensation This might happen if 

there exists space like zero energy excitation of pionic 
modes. The short-range correlation, on the other hand, 
removes such a possibility at least at densities near 
the saturation densities. In the context of relativistic 
heavy ion collisions (RHIC), the importance of medium 
modified pion spectrum was discussed by Mishustin, 
where it was shown that due to the lowering of energy, 
pion, in nuclear matter, might carry a bulk amount of 
entropy . Subsequently, Gyulassy and Greincr studied 
ionic instability in great detail in the context of RHIC 
The production ofpionic modes in nuclear collisions 
was also discussed in 

In experiments medium dependent pion dispersion 
relation can also be probed via the measurements of 
dilepton invariant mass spectrum. The lepton pairs 
produced with invariant mass near the p pole are 
sensitive to the slope of the pion dispersion relation 
in matter Q. Particularly the softening of momentum 
dependence of the pion dispersion relation in matter 
leads to higher yield of dileptons. Gale and Kapusta 
were first to realize that the in-medium pion dynamics 
can be studied by measuring lepton pair productions 
0. Most of the earlier studies of in- medium pion 
properties were performed in the non-relativistic frame 
work [ A [a , |9|. A quasi-relativistic approach was taken 
in [l3, llil, where the calculations were extended 
to finite temperature. In particular, [T^ discusses 
various non-collective modes with the possibility of pion 



condensation. In on the other hand, the dilepton 
production rates were calculated using non-relativistic 
pion dispersion relations. Ref.p^ treated the problem 
relativistically but free Fermi gas model was used, while 
in [T3 | pion propagation was studied by extending the 
Walecka model ll5j including delta baryon. In recent 
years, there has been significant progress to calculate 
dilepton production rates involving pionic pro perties in 
a more realistic framework 0, [12, IH, • 

In the present paper we study pion dispersion relations 
in ANM using relativistic models. This is important as 
most of the calculations, as mentioned above, are either 
restricted to SNM or performed in the non-relativistic 
framework. Here we focus on the propagating modes of 
various charged states of pions which are non degenerate 
in ANM. The importance of relativistic corrections and 
density dependent pion mass splitting in ANM in the 
context of deriving pion-nucleus optical potential was 
discussed in [l^l . The formalism adopted in [l^ was that 
of chiral perturbation theory. Recently, in the context 
of astrophysics, pionic properties in ANM has also 
been studied by invol ving Nambu-Jona-Lasinio model 
[20l . [2H . Motivated by |l9l | present authors revisited the 
problem in ref. [l^l where not only the static self-energy 
responsible for the mass splitting but the full dispersion 
relations for the various charged states of pions were 
calculated after performing relevant density expansion in 
terms of the Fermi momentum. However, in our previous 
work (2^, pions were included via straight forward PV 
coupling in the Walecka model which renders the 
theory non-renormalizable. Although the problem of 
non-renormalizibility could be avoided by considering 
the PS ttN coupling. This, on the other hand, fails to 
account for the pion-nucleon phenomenology. 

Historically, the extension of the Walecka model to 
include the isovector tt and p meson for the realistic de- 
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scription of dense nuclear matter (DNM) while retaining 
the renormalizibility of the theory was first made by 
Serot However, in this work, the calculation was 

restricted only to the mean field level which gives to rise 
tachyonic mode for pions even at density as low as O.lpoi 
where po denotes normal nuclear matter (NNM) density 
[23 |. Such a non-propagating mode for the pions can be 
removed by extending the calculation beyond the mean 
field level as showed by Kapusta 2j|. This, in effect, 
means inclusion of the tt-NN loop while calculating the 
in-medium dressed propagator for the pion. This model 
has an added advantage because of the presence of tt-ct 
coupling in addition to the usual PS coupling of the pion 
with the nucleons which is responsible for the generation 
of small s-wave pion nucleon interaction in vacuum. 
This is consistent with the observed characteristics of 
the pion-nucleon interaction which is dominated by the 
p-wave scattering while the s-wavc scattering length is 
almost zero. In matter, however, as argued in [U, [2^, 
such subtle cancellation does not occur resulting in a 
unrealistically large mass for the pions in matter. To 
circumvent this problem it was suggested in [^J] to use 
the pseudovcctor coupling even though it makes the 
theory non-renormalizable. 

The theoretical challenge, therefore, is to construct 
a model with ttN PV interaction which preserves the 
renormalizibility of the theory. This was accomplished in 
ref. 1251 following the technique developed by Weinberg 
[26l. l27l. [28t and Schwinger [23 • Here one starts with the 
PS coupling and subsequently invokes non-linear field 
transformations to obtain PV representation. Unlike 
straight forward inclusion of PV interaction in this 
method one requires only finite number of counter terms 
which makes the theory renormalizablc. We, here, start 
with this model developed by Matsui and Serot [2^ to 
study the pion propagation in ANM. Clearly, the model 
adopted here is different from what we had invoked in 
our previous work Furthermore, in [22|, for the 

determination of pion self-energy in matter only the 
scattering from the Fermi sphere was considered and the 
vacuum part was completely ignored. The latter gives 
rise to a large contribution to the pion self-energy in 
presence of strong scalar density (ps)- 



The above mentioned model has various shortcomings 
too. In fact; the ref. [1^ itself discusses its limitations 
in describing many body nN dynamics. For example, 
the successful description of the saturation properties 
of nuclear matter in this scheme requires higher scalar 
mass which gives rise to larger in-medium nucleon mass 
compared to the MFT. In addition, it also fails to 
account for the observed pion-nucleus scattering length 
at finite density psj . In the same work, chiral tt-ct 
model has also been discussed to which wc shall come 
later. In the end, wc present results calculated using this 
non-chiral model together with what we obtain from a 
chirally invariant Lagrangian. 

In [12] we have discussed another interesting possibil- 
ity of the density driven n-rj mixing in ANM. However, 
quantitatively, the mixing was found to be a higher order 
effect and does not affect the pion dispersion relations at 
the leading order in density. Hence in the present paper 
we neglect tt-j] mixing. 

The plan of the paper is as follows. In section II we 
present the formalism, where we start with PS coupling 
in subsection A. In II B we invoke non-linear field trans- 
formation [25! and subsequently report results involving 
PV coupling. In section III we present results using 
recently developed chiral effective model in the context 
of nuclear many body problem [sol . [3l| . Finally, section 
IV presents the summary and conclusion. Detailed 
derivations for the Dirac part of the pion self-energy for 
PS and PV couplings have been relegated to appendix 
A and B respectively. 



II. FORMALISM 

A. Model with pseudoscalar ttN interaction 

Wc start with the following interaction Lagrangian 
jiven by [1^, 



I 

where. Here, 4', <i>7r, ^s, '^p and represents the nucleon, 

pion, sigma, rho and omega fields respectively and their 

Gf,^ = - d^'P^p (2a) 

Bp.u = dp,^p^ - d^^pp - gpipp X (2b) 
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masses are denoted by M, m^, rus, mp and rrii^. This 
model successfully reproduces the saturation properties 
of nuclear matter and yields accurate results for closed 
shell nuclei in the Dirac-Hartree approximation [s^ . 

It is to be noted that in Eq.([T]) the pion-nuclcon dy- 
namics is described by 



PS 



1- 



(3) 



where, is the pion-nuclcon coupling constant with 

Ij- — 12.6 [11]. This apart, the interaction Lagrangian 
of Eq. ([1]) also has another term involving the coupling of 
pions with the scalar meson given by 



[3J| which does not distinguish between n and p. Here, 
for the moment we neglect explicit symmetry breaking 
{n-p mass difference) i.e. M* = M* = M*. 

It is to be noted that in the mean field theory (MFT). 
only Fig.([T^), i.e. the tadpole diagram contributes, while 
Fig.lUb) is neglected. The origin of tachyonic mode can 
now easily be understood. The pion mass in matter due 
to the tadpole is given by [l^l 



2 94>^9s 



{Pn + Pp) 



(7) 



(4) 



Here, is the coupling constant of the scalar to pion 
field. The ttN scattering amplitude would now involve 
both nucleon and sigma meson in the intermediate state 
causing sensitive cancellation between the two that gives 
reasonable value of the s-wave scattering length [2J1 as 
mentioned before. At the self-energy level Eq.Q and 
([4]) will generate the exchange and the tadpole diagram 
as shown in Figdb and[T^. 



The second term of the last equation is quite large 
even at densities far below po density viz. for p ~ 0.1 po 
m*J < 0. 
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FIG. 2: (a) represents the one-loop self-energy diagram for 
tt", and (b) represents the same for tt*. 
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FIG. 1: Tadpole contribution to the pion self-energy 



First we consider the tadpole diagram whose contri- 
bution to the self-energy is given by S^'-^ = — 50ir'7is(/>o 
where, <t>Q = and p''{= p^ + pfj. Here pf (i = p, n) 

represents scalar density given by 



M* 
27r2 



E*h - M*^ In 



E* + h 
M* 



(5) 



The effective nucleon mass M* as appears in Eq.® 
can be determined from the following self-consistent con- 
dition [H. 



mi ^ 



(6) 



It is clear from Eq.® that AAf* = M„ - Mp = AM 
as the nucleon masses are modified by scalar mean field 



FigQb would involve various combinations of n and p 
depending upon the various charged states of pions as 
shown in Fig. (Hi) and Fig.fgb). 



S*(g) 



X Tr[{ir{q)}iG,{k + q){ir{-q)}iG,{k)] 



(8) 



where the subscript denotes cither p (proton) or 

n (neutron). F(g) is the vertex factor. F = —175 or 

~*m^'^5'^A'9^ ^'^^ ^^"^ coupling respectively. Ex- 
plicitly, 



G,(fc) = Gf(fc) + Gf(fc). 



(9) 



where, 



Gf(fc) 
Gf(fc) 



M* 



P - M*2 + i( 
E* 



(10a) 

Siko~E*)eikf ~\k\){10h) 



Here, Gf{k) and Gf{k) represent the free and the 
density dependent part of the propagator. In Eq. pITl) 
k is the nucleon momentum; kf denotes the Fermi 
momentum and M* is the in-medium nucleon mass 
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modified due to scalar mean field [3J]. We, from now 
onward, use kp and fc„ to denote the proton and neutron 
Fermi momentum respectively.The nucleon energy is 
E* = yMpTk2. 

Note that the total self-energy is given by ^^gtaiil) — 
+ S'^'^. Using Eq.® and Eq.^, the expression 
for self-energy given in Eq. ([8]) takes the following form : 



Thus, the pion self-energy can now be written as: 



(9) 



^*{FD+DF) 



*DD 



(q) (11) 



Here g is (/tt/wtt) for PS* (PV) coupling. For tt* the 
coupling constant g-^r (or /t) gets replaced by v^^tt- The 
values of the coupling constants and /jr are determined 
experimentally from ttN and NN scattering data. T is the 
trace factor which consists of four parts : 



^FF 



-iFD 



-iDF 



(12) 



Detailed expressions for T^^, T^^ and T^^ will 
be discussed later. Here the term contains the 

product of two delta functions {G^{k)G^{k + q)) which 
put both the loop-nucleons on shell implying the cut 
in the loop (Fig[3^,). This means that pion can decay 
into nucleon- antinucleon (Fig|5}3) pair which happens 
only in the high momentum limit i.e q > 2kp_„ and also 
qo > 2Ep^ where is the Fermi energy for proton (or 
neutron). Under this conditions only contributes 
to the self-energy. But in the present calculation, we in- 
vestigate low momentum (of pion) collective excitations 
only Therefore T^^ {i.e S*^^(g)) is neglected. 



<t>-- 



N 




N 



FIG. 3: (a) represents cutting of loop implied by the product 
of two delta functions and (b) represents the decay of pion 
into nucleon- antinucleon. 



The self-energies for different charged states of pion 
arc calculated using the one-loop diagram shown in the 
Fig. ([2)3). The first term of Eq. (fT3)) is same as the pion 
self-energy in vacuum with Mi M* . This part is di- 
vergent. 



-,FF 

-PS 



2Tr [75zG^(fc)75zG^(fc + (7)] 
k-{k + q) 



M 



(/c2 - M*2) ((fc g)2 - M*2) 



(14) 



Here factor 2 that appears in Tp;^ , i. e. in Eq. (fT4|) , fol- 
lows from isospin symmetry for M„ = Mp. The FF part 
of self-energy for PS coupling is calculated from Eq.lTF 
by substituting Tp^ and it is denoted by H*/'/^ {q). 



d'^k 
M*2 - 



k-{k + q) 



(fc2-M*2)((fc-hg)2-M*2) 



(15) 



From Eq. (fT5|) it is observed that 'S'^g^{q) is quadrati- 
cally divergent. To eliminate these divergences we need 
to renormalize Sp^^('7)- Here we adopt the dimen- 
sional regularization [Sa, [13, [sl] technique to regularize 
Y,*jfg^{q) with the following results (details are discussed 
in appendix A). 
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mi J 





8AP{M - M*)^ 



tan 



q 



r*2 _ g2 



2A/2 /4M2 - m2 
tan 



(4Af2 - ml) 



.2 , ml{M-M*f M , , 

(4M2-m2) mr " 



+ (A/2 - M*^) + 



8Af2 



tan 



^AAP-mi^ 
m^ 



^4A/2 _ 



fix 3x 1 - x)q In —2 2— 

\Ai^ — m^xlyl — X) 



(16) 



It is found that the result given in Eq. (jT6|l is finite 
and no divergences appear further. In the appropriate 
kinematic regime it might generate imaginary part: 



The {FD + DF) part of the self-energy for 7r° and i:^ 
can be written as 



lTni:%^{q) = dx {M*^-'iq^x{\^x)) 

'^^ Jo 

X Im[ln(7\/*2-q2a;(i_2;)_i^)] 

r,2 r 1 

J^/q^~^AM*^ 9 (q^ ^ 4M*'^) 



(17) 

If we consider that (Af * — M) is small enough then the 
term ln[(A/*2_q2^(-;^_^)>)y(-^2_^22,(x_-a;))] ofEq.^ 
can be approximated to 21n(A/*/Af) and the last term 
of Eq. pB)) can be easily evaluated to give 



where, 



^*p'i;{q,m^)^-C+Vq'. 



(18) 



^*±{FD+DF) 
^PS 



where, 



Bps 



-89: 



(Pk 



^PS 



(22) 



[Ap5 =F Bps] 

*0(FD+DF) , j-^*{FD+DF) 



^PS 



>iq) T <5S;^s 



(k ■ qf 



-A{k-qf 
q\k ■ q) 
q^ - 4(fc • g)2 



(23) 



(24) 
(25) 



^ — 27r2 



3(2Af2 - Af*2) + 2Af*2 In 



The trace of (FD+DF) part for 7r°, 



(19) 



+ 75G^(fc + <7)75G^ (fc)] + b ^ 420) 



and for 7r+' \ 



with 



d{k 



■p,n 



|k|). We restrict ourselves in the 
long wavelength limit i.e. when the pion momentum (q) 
is small compared to the Fermi momentum (fcp,,i) of the 
system where the many body effects manifest strongly. 
In this case particle propagation can be understood in 
terms of collective excitation [35| of the system which 
permits analytical solutions of the dispersion relations 
[35I . [3^ . But in the short wavelength limit i.e. when 
the pion momentum (q) is much larger than the Fermi 
momentum (fcp^„), particle dispersion approaches to that 
of the free propagation. Note that for SNM Bps = 

1 • ^*±(FD+DF) ^*0(FD+DF) 

implymg Epg^ = Sps 



rpFD I rpDF rp 

Tps +Tps = Tr 



75Gp(„)(fc + (?)75G,f(p)(fc) 
+ 75G^(„)(fc + g)75G^(p)(fc)j (21) 



In the long wavelength limit we neglect the term q"^ 
compared to the term 4(fc • q)"^ from the denominator of 
both Aps and Bps in Eqs.(l24|) and ([251). Explicitly, 
after a straight forward calculation we get, 
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^*0{FD+DF) 
■"PS 





)1 







(26) 



and 



c^*{FD+DF), ^ 



27r2 
27r2 



Co + Vp 



1 



Co 
Co 



M* 



tan 



tan 




(27) 



where i;p,„ = kp^n/E* „, E* „ = ^ M*^ + fc^ „ and cq = 

qo/|q|- The approximate results of Eqs. ipS]) and ([271) are 
given below. 



-A-B-T+g 



(28) 
(29) 



where, 



*2 



*l ~ [nps + tF"^ + mlo] 

and 

Ops + 



1 T 5^ps 



(34) 



where, 



B 

f 

Q 
S 



a. 



1 (K 

3 \ B* 



M* 



SI'S 



9t 



1 fii 

5 \e;^ 

1 

3 \^M*2 



*2 



A/* 



The self-energy for PS coupling: 



E^°/(g) = Sj4(,,m.) + Y:^,-'-^->{q) (31) 

The dispersion relations can be found by solving 
Dyson-Schwingcr equation. 



*0,±(FD+DF) I 



(30) 



nps ^g-A-B-T 

5^PS = 



(35) 



Now we presents the dispersion relations for 7r°'^ with 
the effect of Dirac sea. 



(36) 



The effective masses (m*?,±) with Dirac sea for differ- 
ent charged states of pion are given by 



where. 



*2 



*2 



{Aps-mlo)/V 
and 



(A 



PS -m 



(1 T <5Aps)I? 



(37) 



m^o,±-(E*"'±(g) + S^^) = 



(32) 



Here m^o,± are the masses of 7r° and tt^. The disper- 
sion relations without the effect of Dirac sea for tt'^'^: 



The effective masses without Dirac sea are 



(33) 



Aps =C- Qps - 
SAps = 



(38) 



The PS coupling the asymmetry driven mass splitting 
is of 0(fc3(„)/M*2). The terms SAps and S^ps are 
non- vanishing in ANM and responsible for the pion mass 
splitting. 
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TABLE I: This table presents the effective pion masses includ- 
ing the tadpole contribution to the self-energy in PS coupling. 
Kapusta corresponds to ref. [2^ and BDM corresponds to the 
present calculation. 



m*3) 



*2 



MFT mlo + E^"^ 

Kapusta ^ps + (S^'^ -|- m^o) 
BDM [C-{9.ps + Yi^'^ +mlo)]/T) 



In Fig|3]and[5]wc present the density (p) and asymme- 
try parameter (a) dependent effective masses for the var- 
ious charged states of pion. In the top panel we present 
the results without vacuum correction (Dirac sea). Here 
we include both the tadpole and n-n loop. 



800 
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1.0 1.2 1.4/ 1.6 1.8 2.0 

P/Po 



FIG. 4: Nuclear density (p) dependent effective pion masses 
for PS coupling at a = 0.2. The dotted, dashed and solid 
curves representing tt", tt"*" and tt'' without (upper panel) 
and with (lower panel) the Dirac sea effect. 



It is evident that the inclusion of ^p) removes the 
tachyonic mode but gives rise to effective pion masses 
which are unrealistically large as discussed by Kaputa 
[2^ as shown in the top panel of FigH) 

It is to be noted that the inclusion of the vacuum part 
reduces the effective pion masses and gives reasonable 
value for the density dependent pion masses in matter at 
NNM density. The reason for this could be understood 
from the TableUwhich enumerates expressions for the ef- 




150 



0.0 0.2 0.4 0.6 0.8 1.0 

a 



FIG. 5: Asymmetry parameter (a) dependent effective pion 
masses for tt" (solid curve), tt"*" (dashed curve) and 7r~ (dotted 
curve) in ANM at NNM density for PS coupling. The upper 
and lower panel represents effective pion masses without and 
with vacuum correction. 



fective pion masses that we obtain in three different cases. 
The top row represents effective pion masses for the case 
considered in [23[ which gives rise to the tachyonic mode, 
the second row corresponds to the case discussed by Ka- 
pusta [iSl and in the last row we present results of the 
present work as by BDM. The presence of the additional 
term 2? somewhat tames the dispersion curve bringing 
the masses down compared to j24| . This can be noted 
that at the MFT level S-^*^ involves sum of the scalar 
densities and Pp. Therefore, in MFT, as expected, 
the masses are insensitive to asymmetry parameter a. 



For Pb-like nuclei (a = 
Am^+ = 17.37Afel/ and 
vacuum correction. 



0.2), Am^o = 16.8MeV, 
Am^- = llAlMeV with 



The dispersion relation for various charged states 
of pion are shown in the FiglS] where upper and lower 
panel present the dispersion curves without and with 
the Dirac sea including the tadpole contribution. In the 
presence of Dirac sea 7r+ and tt^ dispersion curves are 
indistinguishable. 
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nonlinear chiral transformation |2£ 



^ = 



1 - i75T • ^ 



i+e 



(39a) 



M ~ g.'f, + {M - gs'fsr + gin 



(39b) 

(39c) 



The last two equations (|39bp and (|39cp are used to 
express the old fields (f>s and in terms of new fields 
and $' . 



FIG. 6: The dispersion relations of tt" (sohd curve), 7r+ 
(dashed curve) and tt" (dotted curve) for PS coupling at 
p = 0.17/m~^ and a = 0.2. The upper and lower panel 
representing pion dispersions without and with the Dirac sea. 



B. Renormalizible model with pseudovector nN 
coupling 

To obtain the pseudovector representation from the 
starting Lagrangian given by Eq.((T]), one gives following 

I 

The transformed Lagrangian is 

C = i''(i7,,a^-M)*'-ig^^''7^(f-$^)*' + .g,,*'$;*'-5^§'7,,$^*' 
+ i(a^$,a^$, - ml^l) + i(9^$, - .gp<gp^ X $,) • (a^<l, - 5p$^ X 

^2 



1 - 2(/Jm,)$; 



1 + (/./m.)2$;2 
(1 - {U/m^r$'^^)^i + (s.M)(/./mJ$;2 



(40a) 



1 + (A/m,)2$;2 



(40b) 



(/tt/TOtt)^ 



1 + iu/m^yK' 



*'7p(r-^;)x(a^$; 



1 + {U/m^r<^'^ 



*'757MT-(a^^;-.9p^;:x$;)vi/' 

(41) 



We see from the Eq. ST] that the n-N PS coupling has 
disappeared and instead the pion-nucleon dynamics is 
now governed by the last term of the above mentioned 
equation. At the leading order one obtains the usual PV 
coupling represented by, 



Z:^^ = -^*'757^5p (f . $;) ^' (42) 

Here /tt is the pseudovector coupling constant and 

|j = 0.08 [i^l. First we discuss the FF part where the 
trace factor is given by 
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rf^FF 



-2 Tr[^5^''q^iG^{k)^^-rq,iG^{k + q)] 

'M*^q^ + k-{k + q)q^ - 2{k ■ q){k + q) ■ q 



(fc2-M*2)((/e + g)2-M*2) 



(43) 



Now the FF part of the sclf-cncrgy for PV coupling is 
denoted by Spv^(g). From Eq.^ and Eq.^^ we get, 



8i 



d'^k 



m^j J (27r)4 
M*2g2 + A: . (/fc + q)q^ _ 2(fc • q){k + q) ■ q 



{k"^ ~ M*^){{k + qf - M*^) 



(44) 

Direct power counting shows that the term T,*Jlf'{q) 
is divergent. The appropriate renormahzation scheme 
for the present model has been developed in ref[25j. We 
first consider a simple subtraction scheme described in 
Appendix B to obtain 



i*_R / \ 9^ f fir 



2M 



*2 



dx In 



/ A/*2 - q^x{l - x) 
VAf*2 - m2x-(l - a-) 



Now Epy(g) can be approximated to 



i:fy{q)c^Cq'-Vq'' 



(45) 



(46) 



where, 



V 







2 










f 











(47) 



The FF part can also develop imaginary part as given 



by 



im (q) 



2Af*V9^ - 4A/*2 e {q^ - iM*^) 

(50) 



It is observed from Eq.dSQl) that Im T,*JlF{q) IS non- 
vanishing only if g2 > 4A/*2. 

The trace of the (FD+DF) part for 7r°, 



rpFD , rpDF 

^PV + ^PV 



and for 7r+'^ \ 



Tr[-/5^G^{k + q)j,^G^ik) 
l5iG^{k + qh5iG^{k)]+[p^n] 



(51) 



FD 
PV 



DF 
PV 



Tr 



75^G^(„)(fc + g)75^G^(p)(fc) 



+ ldG^in){k + q)idGl^p){k)\ (52) 

In pure neutron (or proton) matter one of the terms of 
Ea. ((52l) viz G^^^j = for the charged pion states. The 
same argument holds true for the neutral pion where only 
two terms would contribute which can be observed from 
Ea. (j5ip . In case of pure neutron (or proton) matter p{n) 
appears as the intermediate state. Now the [FD + DF) 
part of the self-energy for tt" and tt^ can be written as 



^*0{FD+DF) 
^PV 



^*±{FD+DF) 
-'PV 



(9) = -8 ( ^ 



d^k 



{2ttYE' 

d^k 
{2tt)^E' 



i-PV 



(53) 

■[Apv =F Bpy] 



,~^*OiFD+DF), ^ c^*{FD+DF) f s 



(54) 



On the other hand borrowing results from [23 one has. 



where 



where. 



C 
V 



^*P^y{q)c^C' +V'q^ 

{lt^)\mM - M*)ml] ^ 
-h\\2M*Hn{M*/M)\ _ 



(48) 



(49) 



It might be mentioned, although C, I? arc different from 
C'jV, their effect on the effective pion masses and cor- 
responding dispersion relations are found to be marginal 
as we discuss later. 



Apv 



M*'^q^ 



q^ - 4(fc • g)2 

4M*2g2 



1 



q^ - 4{k ■ qf- 



p + 0n) (55) 

(fc-Q)(0p-0„) (56) 



In the long wavelength limit considering collective ex- 
citations near the Fermi surface, {FD + DF) part of the 
pion-self energy can be evaluated analytically. In this 
case wc can neglect the term q^ compared to the term 
4(A: • (7)2 from the denominator of Apy and Bpy in Eqs. 
([55|) and (|56p . This is called hard nucleon loop (HNL) 
approximation [s^ . Explicitly, after a straight forward 
calculation, we get. 
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and 



In 



Co In 



Co + Vp 



Co - Vp 



In 



l-v„ 



Coin 



Co +v„ 



Co - Vn 



(57) 



ST. 



*{FD+DF) 
PV 



2 , , A'/*2g2 / 



The approximate results of Eas. (|57|) and ([55]) are 



■ (FD+DF) 
PV 



{q) ~ £qo 



(59) 
(60) 



where, 



B 



711-!^ IT 



> 


( kl 












( fcp 




> 




+ *) 




(A 




j 







(61) 



The total pion self-energy for PV coupling is 



Co + Vp 



Co - Vp 

Co + Vn 



2M* 



Co - Vy, 



2M* 



tan" 



-1 ^py^ 



tan 



Co A/* y 

-1 fcn'v/Co — 1 



coM* 



r 



1 



f2f 



(58) 



(65) 



These results are the same as that of [2,21 with some no- 
tational difference such as fipy fi^^, SVLpy SVt^^, 
fipv/(l - Vtpv) Xtttt and B/{1 - fipy) Pt,-^. In 
([57)1 and (5Apy and 6V,py are responsible for the 

asymmetry parameter (a ~ ''p" +p'' ) dependent mass 
splitting, where and Pp are the neutron and proton 
density respectively. Clearly for SNM 6Apv and 6ilpv 
vanish. 



^*0,±/ N / \ , ^*0,±(FD+DF) , . 

Spy [q) = ^PVW + ^PV il) 



(62) 



The approximate dispersion relations and the effective 
pion masses of different charged states in ANM without 
and with the Dirac sea effect are presented below. 

The dispersion relations for tt^'^ without the effect of 
Dirac sea are as follows. 



7^ 

ITTTT 



m*2 



(63) 



where '7i*o,± is the effective pion masses without Dirac 
sea effect: 



The dispersion relations for tt'^'^ including the effect 
of Dirac sea are given by. 



r^2 



277r7r) 



(66) 



The effective masses (m*) of different charged states 
of pion are found from Eq. (j66|) in the limit |q| = 0. 



*2 



1 - n 



and TO* ± 



9 

TO^ + 



PV 



1 - {Qpv±Snpv) 



(64) 



*2 



1 - A 



and mil 



PV 



1 - {Apv ± 6Apv) 



(67) 



where. 



where. 
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A 



PV 



A + B + C 



SA 



PV 



1 - 



1-Ap 



1-Af 



V 
1-Ap 



(68) 



This is to be noted that, if one use Eq.(|^51) instead of 
Ea. (|46)) : Ti^o.i and C will be replaced by (to^o,± + C) 
and V respectively. (Jtttt will vanish. Numcrically,as 
mentioned before. Eg. pS)) and Eq. (|15)) give results 
very close to each other. Clearly from Eq. (|CT|) . £ 
indicates that the asymmetry driven mass splitting is of 
0(fc5^„)/M*4) for PV coupling. 
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FIG. 7: Effective masses (for PV coupling) of tt" (solid curve), 
tt"*" (dashed curve) and tt~ (dotted curve) are represented 
without Dirac sea (upper panel) and with Dirac sea (lower 
panel) at a = 0.2. 



To quote typical values of the pion mass shifts for PV 
coupling at normal nuclear density (po = 0.17 fm^^) 
for Pb-like nuclei which are Am^o — 6.07MeV, 
Am^+ — A.6MeV and Am^- = 8.Q2AIeV with vacuum 
correction and the corresponding values are A.QbMeV, 
SAlMeV and 6.82MeV without vacuum correction. 

In Fig[7|we show results for the density dependence of 
effective pion masses for various charge states at a = 0.2. 



It is observed that the 7r~ mass increases in matter 
while 7r+ decreases at higher density. The mass splitting 
is quite significant even at density p > 1.5po- In the 
lower panel we present results with vacuum corrections. 
Evidently the effect of vacuum corrections is found to be 
small. 
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FIG. 8: Asymmetry parameter (a) dependent effective masses 
(for PV coupling) of tt*^ (solid curve), tt^ (dashed curve) and 
7r~ (dotted curve) at p = 0.17/m~^ without (upper panel) 
and with (lower panel) vacuum correction. 



It should however be mentioned that the vacuum 
correction part for PV coupling is rather small. For 
loops involving heavy baryons it could be quite high. 
For detailed discussion we refer the readers to [U, \^ . 
In present case we take only the nucleon loop in presence 
of the scalar mean field. 

We also present results of asymmetry parameter de- 
pendence effective masses for different charge states of 
pion in Figl8]at normal nuclear matter density. The up- 
per and lower panel present the effective pion masses 
without and with vacuum correction. It can be ob- 
served that the asymmetry parameter dependent pion 
mass splitting is insensitive to the vacuum correction. 
The pion dispersions in medium for various charge states 
of pion are presented in Figl9]for PV coupling. 



III. MODERN TECHNIQUE 

In the previous sections we have discussed pion prop- 
agation in ANM using both the PS and PV interaction 
within the framework of non-chiral model. However, 
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FIG. 9: Pion dispersion relation without (upper panel) and 
with (lower panel) the effect of Dirac sea for PV coupling. 
The solid, dashed and dotted curves respectively indicate the 
dispersion curves of tt", tt"*" and 7v~ at p = 0.17/m~^ and 
a = 0.2. 



mentioned before, it fails to account for the pion-nucleus 
dynamics in nuclear matter both in the PS and PV 
representations. In fact, it gives too strong pion nucleon 
interaction in matter which cannot be adjusted by fixing 
the s-wave tt-N interaction in free space even in PV case. 
In this context the Dirac vacuum involving baryon loops 
was found to play a significant role. If one uses the chiral 
model and breaks the symmetry explicitly, the results 
are found to be very sensitive to the renormalization 
scheme [l^. In [i^ it was shown that the relativistic 
chiral models with a light scalar meson appear to provide 
an economical marriage of successful relativistic MFT 
and chiral symmetry. It, however, fails to reproduce 
observed properties of finite nuclei, such as spin-orbit 
splittings, shell structure, charge densities and surface 
energies. Since then, there has been series of attempts 
to construct a model which has the virtue of describing 
both the properties of nuclear matter and finite nuclei 
[13, HH, Sll, m, m, Si], currently, the non- linear chiral 
effective field theoretic approach seems to be quite 
successful in this respect. It might be recalled here, that, 
in such a framework, the explicit calculation of the Dirac 
vacuum is not required, rather, on the contrary, here, the 
short distance dynamics are absorbed into the parame- 
ters of the theor y adiust ed phenomenologically by fitting 
empirical data [30, [3l|, S^l- ^ow we proceed to cal- 
culate the effective pion masses in ANM in this approach. 



the interactions as represented by Eq[l] and Eql41] fail 
to describe in-medium t:N dynamics as shown in [25l |. 
It was also observed that the chirally symmetric model 
(linear) has also various limitations [25l |. For example, as 



By retaining only the lowest order terms in the pion 
fields, one obtains the following Lagrangian from the chi- 
rally invariant Lagrangian[3lj : 



Jtt ^ 



1 



NL 



(69) 



The terms Cnl and SC contain, respectively the 
nonlinear terms of the meson sector and all of the 
counterterms. The explicit expressions for Cnl and SC 
can be found in 31 1 . 



It is to be noted that the meson self-energy can be 
found by differentiating the energy density [3l| at the 
two-loop level with respect to the meson propagator 
as indicated in Fig (TO] One may therefore, identify 
the FF, FD and DD parts of the self-energy with the 
vacuum-fiuctuation(VF), Lamb-shift(LS) and exchange 
(EX) contributions to the self-energy respectively. The 
VF and LS terms are related to the short-range physics 
while EX part is related to the long-range physics. The 
detailed discussion about this short and long distance 
separation can be found in [s^, HH, H^]. The diverging 



FF part of the self-energy and LS can be expressed as 
a sum of terms which already exists in the effective field 
theoretical Lagrangian and can be absorbed into the 
counter terms. The short distance physics, as shown 
in [3l|, while calculating exchange energies, are either 
removed by field redefinitions or the coefficients are 
determined by fitting with the empirical data. The 
long-range part is computed explicitly that produce 
modest corrections to the nuclear binding energy curve. 
This can be compensated by a small adjustment of the 
coupling parameters. 

Recently in ref.(3l| the exchange energy contributions 
of pion has been calculated within this theoretical frame- 
work. We adopt the same parameter set as designated 
by MOA in [311 to calculate the tt self-energy explicitly. 
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FIG. 10: Diagrams correspond the two loop self-energy and 
energy 
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FIG. 11: Effective pion mass at different densities with a = 
0.2. 



The corresponding results are presented in Fig llll Here 
we simply depict the final results as the expressions, at 
this order, for the pion self-energy and density depen- 
dent masses of 7r° and tt^ remain same as those of Eg 1571 
except for the coupling parameters. Quantitatively, it 
is found that, for the lower density, i.e p ~ poj the ef- 
fective masses for 7r~, tt*' and t:~^ states are comparable 
with that of PV coupling (FigO, while at higher density 
the mass splitting is significantly enhanced. The charged 
states, i.e. tt^ show stronger density dependence com- 

I 



Here e ~ 2 — ^■ and p, is an arbitrary scaling 
parameter, is the Euler-Mascheroni constant and 



pared to PV coupling. We also observe that the density 
dependence of tt*' is rather weak. 



IV. SUMMARY AND CONCLUSION 

In the present paper pion propagation in ANM 
has been studied within the framework of relativistic 
hadrodynamics in presence of the scalar mean field. 
We start with the model developed in (2^ and present 
analytical results for the pion dispersion relations in 
ANM by making HNL approximation and suitable 
density expansion. Su bseq uently, we invoke the chirally 
invariant Lagrangian [47l . l48l | by retaining only the 
lowest order terms in pion field and compare the results 
with non-chiral model calculations performed in section 
II. 

The splitting of the various charged states of pion even 
at normal nuclear matter density is found to be quite 
significant. Such mode splittings in ANM is, infact, a 
generic feature of all the isovector mesons. Therefore, 
it would be interesting to estimate similar splitting for 
the p meson and other isovector states. It is to be noted 
that the mass s plit t ing is related to the pion-nucleus 
optical potential [1^, [23 . As for the dispersion relations, 
we restricted our calculation only to the time-like region 
which can also be extended to study the space like 
modes both for the pions and rho mesons. 

Acknowledgment : The authors gratefully thank 
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suggestions. 



APPENDIX A 

After using Feynman parametrization, the term 
T,*pg^{q) in Eq. p5|) can be written as 



(Al) 



7^ = (7^; — 1). The imaginary part of Y,*^^ {q) can easily 



^*FF 
^PS 



(9) 



27r2 



9i r 



JO 



dx (Att^ 



27r2 



#i / dx(M 



A. 1 

27r2 e 

*2 



dx 

m 

1 - e 

il/*2 - 



M* 



k-{k + q) 



{{k + qxY + q^x{l - x) - M*^Y 
"M*2 - 3g2x(l - x) + 2eq^x{l - 



(A/*2 - q^x{l - x)y 



9l 



■iq^x{l 



27r2 
)) In [M' 



M 



*2 



(7^-ln(4V)) 



q^x{l-x)) 



V 
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be found by simply replacing In (M*^ — q'^x{l — x)) 
with In (M*^ — q'^x{l — x) — irj^ where 77 is an arbitrar- 
ily small parameter and the term iij comes from the 
denominator of Gf when Feynman paramctrization is 
performed considering iC, in the denominator of the 
propagator. 



It is clear from the expression of Eq. (jAip that the sec- 
ond term is divergent in the limit e — > (as ^ 4). 
To remove the divergences we need to add the countert- 
erms [25| in the original Lagrangian interaction . The 
diverging part of Eq. (jAl[) is 



Here the term In (A/*^ — q^x{\ — x)) has branch cut 
only for M*^ - q^x{l~ x) < and it begins at q"^ = AM*"^ 
i.e. the threshold condition for nucleon-antinucleon pair 
production. So the limit of x-integration changes from 

(0, 1) to (i — ia, i + ia) where a = y^l 
used Imln (Z — irj) = — tt. Now, 



4M' 



and we 



dx e (q^ - 4Af 



*2\ 



AM* 



Now the imaginary part of (q) is, 



(g2 _ 4A/*2) 



(A2) 



V 



PS 



9l 




27r2 




9l 


'AP 2 


27r2 


e e 



^2 j,2 



2e 



(A4) 



l^Y.*/s^{q) = /' dx {M*^-?>q^x{\~x)) 

X Im[ln(7\r2-q2a;(l-x)-i77)] 
^2 



47r 



(A3) 



In Eq. (|A4p we substitute the effective nuclcon mass 
M* = (A/ — g 1(1)0) where M is the nucleon mass and c/jq 
is the vacuum expectation value of the scalar field (j)s. 
The expression given in Eq. (|A4[) tells us that we need 
to be added four counter terms [2^ with the original 
interaction Lagrangian to remove the divergences from 
Sp^^. Therefore the counter term Lagrangian [2^ is 
denoted as 



(A5) 



r 



The value of the counterterms f3i, /?2, (3^ and /?4 are 
determined by imposing the appropriate renormalization 
conditions. 



Pi 
f3i 



\ dq^ 



-5X + g. 



dM 

2 (d^^'p'sil) 



q^—ra 



(A6) 

(A7) 
(A8) 

(A9) 



G-n ~ [<f ^'>n'i^^*psi^)\^^ reproduces the physical mass 
of pions in free space. The counterterm /?4 determines the 
strength of coupling of the vertex. In fact Ep^ (g) 

is found by simply replacing M* with M in Eq. (|Aip . 
We can set (5A = to minimize the effects of many-body 
forces in the nuclear medium ^25i] which is consistent with 
the renormalization scheme for scalar meson. Using the 
conditions given in Eqs. (|A6P - (jA9[) the following results 
are found : 



Here (3i and /32 are the wave function and pion mass 
renormalization counterterms respectively while (3^ and 
/34 are the vertex renormalizaton counterterms for the 
0s $2 vertex and ^^^^ vertex respectively. The condi- 
tions of Eq. (jA6P - (|A7p implies that the pion propagator 
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Pi 



A. 

27r2 
27r2 
27r2 



3 2 V e 



9l ' 



7^ + ln(47r^2 
3a;(l - x) In (M^ - mlx{l - x)) 



dx 



NPx(1 - x) - ■iralx'^{l ~ x) 



\21 



M2 - m2a;(l - x) 



(AlO) 



P2 



2 r 9 
mt 



27r2 



5? ^ 



- 



-7^ + ln(4V) 



27r2 



2 



(M^ ~ 3m2a;(l - x)) 



X ln(A/2 - m^x(l - a;))] 



/?3 



27r2 



g,(2M) --7^ + ln(4V) 



+ 



27r2 
27r2 



5s(2Af) / dxln{M^ - mlx{l - x)) 



/■,,rN / , /M2-3to2x(1-x) 



(All) 



(A12) 



Now the renormalized Yj*j^^{q) is 



/32 - /?30O - 77 



ml) 



Ml 



(A14) 



P4 



(72 o2 

-2^^-''^ + 2^ 



2g,s 



27r2 
27r2 



-7B + ln(4V) 



2c/2 / rfxln(A/2 - m2x(l - x)) 
^0 



29I 



U'Pmlx{l ~ x) 
Jo ""^ (M2 - m2a;(l - x))^ 



dx 



(A13) 



Substitutin g Sp^^ (g) from Eg (jAlj) and ^1, /J2, /^s, 
/?4 from Eqs. (|A10p - (|A13P in Eq. (|A14p it is found that 
divergences in Sp^^ (q) are completely eliminated by the 
counterterms. After simplification I]p^(g, m,r) reduces 
to 



27r2 















m2 y 





2Af*2V4M*2 - g2 

tan 



^4Af*2 - g2 



+ 



2A/2y4M2' 



(4A/2 - ml) 
m-iT 



tan 



,,0 m2(Af-A'/*)2 Af ■ , 



da; 3a;(l — x)q In 



(4A//2 - m2 ) ml 
( M*^ -q^x{l-x) 



8Af2 



-V4M2 



tan 



^4Af2 - ,7l2^ 



^4A/2 - to2 



VA'/2 _ mlx{l - x) 



(A15) 



APPENDIX B 



After Feynman parametrization Eq. (j44p reduces to 
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8^ ( ^ ) fi^' 

m. 



(2^ 



dx 



(M*2 + q^x{l -x) + e) ~ 2{k ■ qf 



2 „1 



71- / JO 



((/fc + qxf + 922;(1 - a;) - M*2)2 



27r2 \ 



,2 / r \ 2 r 



2Ar2 (7B-ln(4V))] 



2M* 



('A 

27r2 \ TO^ 



2A/" / dx \n{M*'^ -q^x{l-x)) 




(Bl) 



The imaginary part of Sp^^ (g) can be found as 



Im 



^2A/*Vg^ -4A/ 



*2 



(g2 _ 4M*2) 
(B2) 



It is clear from Ea. (|B2|) that ImEp^^(g) vanishes for 
q2 < 4M*2. With the same argument as stated for PS 
couphng, we excluded the imaginary part. The diverging 
part of T.%^{q) is 



Here we use simple subtraction method to remove the 
divergence. So, the finite FF part of the self-energy is 



^PV 



id) 



27r2 \ rn^ 



2M 



*2 



' dx In AIl^^Ml^A 

\M*2 — m2a:(l — x) 



(B4) 



PV 



27r2 





■2M*2" 




e 



(B3) 
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